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Topological phenomena in gauge theories have long been recognized as the driving force for chiral 
symmetry breaking and confinement. These phenomena can be conveniently investigated in the 
semi-classical picture, in which the topological charge is entirely carried by ( ant i-) self- dual gauge 
configurations. In such an approach, it has been shown that near the critical temperature, the non- 
zero expectation value of the Polyakov loop (holonomy) triggers the "Higgsing" of the color group, 
generating the splitting of instantons into N c self-dual dyons. A number of lattice simulations 
have provided some evidence for such dyons, and traced their relation with specific observables, 
such as the Dirac eigenvalue spectrum. In this work, we formulate a model, based on one-loop 
partition function and including Coulomb interaction, screening and fermion zee modes. We then 
perform the first numerical Monte Carlo simulations of a statistical ensemble of self-dual dyons, as a 
function of their density, quark mass and the number of flavors. We study different dyonic two-point 
correlation functions and we compute the Dirac spectrum, as a function of the ensemble diluteness 
and the number of quark flavors. 



I. INTRODUCTION 

Topological phenomena in gauge theories 
have been discovered more than three decades 
ago, and remain the subject of intense theo- 
retical research ever since. In particular, mag- 
netic objects (monopoles) have been identified 
as a possible source of confinement [TJ [2] , while 
instantons have been proposed as the driving 
mechanism for chiral symmetry breaking [4j [5] . 

The index theorem establishes a direct con- 
nection between the vacuum topology, and 
zero-eigenvalue solutions of the Dirac equation, 
i.e. the so-called fermionic zero-modes. These 
quark states are insensitive to any perturba- 
tive fluctuation of the gauge field, hence encode 
purely non-perturbative QCD dynamics. Fur- 
thermore, lattice simulations have shown that 
the Dirac eigenstates with near-zero eigenval- 
ues — also known as the "zero- mode zone" 



(ZMZ) — directly correlate with local fluctu- 
ations of the topological charge density. Af- 
ter filtering out quantum fluctuations, lattice 
fields reveal nearly (anti) self-dual smooth fields 
responsible for topology and ZMZ states [3]. 
Using only fermionic states attributed to the 
ZMZ (a tiny subset of Dirac eigenstates, of only 
about ~ 10 -4 of all eigenstates) one finds the 
correct pion mass, quark condensate as well 
as many other hadronic properties. On the 
other hand, filtering out the ZMZ states re- 
moves the chiral symmetry breaking and leads 
to drastic changes in the hadronic spectrum 
computed on the lattice. In particular, some 
masses get shifted by as much as ~ 30% and 
parity-doublets appear (for a recent analysis, 
see e.g. [6]). 

This body of results coherently support a pic- 
ture in which the non-perturbative chiral dy- 
namics in vacuum is mediated by instantons. 
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Indeed, instanton model calculations (for a re- 
view see [7 ) , have been very successful in repro- 
ducing the mass and electro-magnetic structure 
of pions [8 , vector mesons [9] , nucleons [TQHT5] 
and even the AI = 1/2 rule for hyperon [T6] 
and kaon [17] non-leptonic decays. 

In the instanton picture, the width of the 
ZMZ depends on the size of the typical "hop- 
ping" matrix element of the Dirac operator be- 
tween two instantons, which is of the order 
~ p 2 /R 3 ~ 20 MeV, where p is the typical in- 
stanton size and R the typical inter-instanton 
density [4]. This value is comparable to the 
typical light quark masses used in may lattice 
simulations, and this explains why the corre- 
sponding results display significant deviations 
from the naive chiral perturbation theory pre- 
dictions. Furthermore, the specific shape of the 
density of eigenvalues p{\) in the ZMZ depends 
crucially on the theory parameters, such as the 
number of light fermions Nf. 

In this work, we will further investigate topo- 
logical phenomena in the semi-classical pic- 
ture, focusing on temperatures close to those 
at which the expectation value of the Polyakov 
line 

(P(x)> = (exp ^ ij*dx* AZ(x, Xi ) y 

(1) 

drastically changes from 1 to 0. 

The gauge invariant expectation value ([I]) de- 
fines the holonomy of the gauge connection cor- 
responding to a full circle around the periodic 
time direction and is related to the free-energy 
F q of a single static quark: 

(P)~exp(-F g /T). (2) 

Hence, a drastic suppression of (P) reflects the 
onset of the confinement phenomenon and we 
shall denote the corresponding critical temper- 
ature with T c . 

From Eq. it follows that a non-trivial 
holonomy (i.e. (P(x)) ^ 1) reflects a non- 
vanishing vacuum expectation value (VEV) of 
the A 4 (x, x±) component of the gauge field. 
In the semiclassical picture, this condition is 
not fulfilled by the standard instantons and 
calorons (finite temperature instantons). In- 
deed, in such topological classical solutions, the 
A4 component of the gauge field vanishes at 
spatial infinity, limi^^^ Ai(x, X4) = 0. 



Classical gauge configurations, generalizing 
the instanton solution to a non-trivial holon- 
omy case, are called the Kr aan- van- Baal- Lee- 
Lu (KvBLL) calorons [H [19]. These au- 
thors have shown that instantons basically split 
into certain substructures, the N c "constituent 
dyons" . 

The 4N C collective variables for the con- 
stituent dyons are different from that of stan- 
dard instantons: while each instanton is iden- 
tified its position, size and color orientation, a 
dyon is specified only by its position and its 
Abelian phases. The solution and the measure 
for the KvBLL configurations correctly returns 
to the standard instanton one in the limit of 
vanishing holonomy. For review of instanton- 
dyons see [20], [21] and references therein. 

In general case, the holonomy can be 
parametrized in terms the VEV of A4 at spatial 
infinity, according to the usual Cartan subalge- 
bra notation. Such an asymptotic A4 field can 
be taken to be spatially constant and diagonal 
in color space, 

A 4 (oo) = (2ttT) diag(/ii... fi Nc ); ^/ij = 

i 

The latter condition follows from imposing a 
null trace. 

Inserting such a constant field plus and its 
quantum fluctuations into the gauge action, one 
finds that the Yang-Mills commutator gener- 
ates mass terms for all non-diagonal compo- 
nents of the gluon field. On the other hand, di- 
agonal components of the quantum gauge field 
commute with the holonomy and remain mass- 
less. 

It is convenient to adopt a cyclic-symmetric 
notations such that Unc+i = 1 + Mi and intro- 
duce the differences between subsequent holon- 
omy parameters, 

Vi = IH+i ~ Vi\ Vi = 1 - ( 4 ) 

i 

These differences determine the masses of the 
non-diagonal (i + 1, i) components of the gauge 
field. Since the dyon cores are made of such 
field components, the i-th dyon has the core 
size ~ \/(2irTvi). On the other hand, the long- 
range dyon interaction is due to massless diag- 
onal fields and is governed by their electric and 
magnetic charges (see e.g. Ref. [21] ). 
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In this paper, we restrict our attention to the 
SU(2) gauge group. In this case, the holonomy 
parameterization simplifies, since there is only 
one diagonal generator r 3 and a single holon- 
omy parameter, /i = /ii = — /i2- In Fig|l] we 
show the locations of the holonomy eigenvalues 
at two different temperatures. At high T (case 
(a)) the holonomy is close to a zero phase, so 

= 2 ji is small. Respectively, the mass of the 
first dyon, called M, is also small. The mass 
of the second dyon is proportional to the com- 
plementary part of the phase circle ^2, which 
at high T takes nearly the whole circle. The 
corresponding L dyon is thus heavy. 

As T — >• T c , the holonomy moves toward 
±1/4 of the circle, as shown in case (b): and in 
this case both dyons have the same mass. This 
position corresponds to phases Hi = ±1/4 and 
the Polyakov line vanishes P — ^ cos(7r/2) = 0, 
which physically corresponds to the confine- 
ment of heavy charges (quarks). 

In Table I we list the quantum numbers of 
the possible dyon solutions for the SU(2) gauge 
group. In such a gauge theory there are 4 
distinct dyons, denoted with M, M, L and 
Z, which are characterized by different (color- 
) magnetic and electric charges. For the gen- 
eral SU(N C ) case there are Mi, M 2 ...Mn c so- 
called "static" dyons with all diagonal charges 
and just one so-called "twisted" dyon L. 

From the technical point of view, the statis- 
tical physics of the dyon ensemble is quite chal- 
lenging compared to that of the "instanton liq- 
uid" [7]. This is mostly because, on top of topo- 
logical phenomena induced by the zero-modes 
of light fermions, one has to deal with long- 
range Coulomb forces and even linearly growing 
potentials due to a "screening" effect, as will be 
discussed below. 

A first, qualitative discussion of dyonic en- 
semble has been performed by T. Sulejmanpasic 
and one of us [22] . It was pointed out that a dy- 
onic system can realize several distinct phases, 
ranging from a gas of individual LL clusters at 
high T, to a strongly-coupled liquid plasma at 
T « T c . Ref. [22] also includes an extensive dis- 
cussion of the lattice-based phenomenology re- 
lated to self-dual dyons. Extension of that work 
to the case of adjoint fermions will be given in 

The present paper can be considered as a 
quantitative continuation of that study, we shall 



T>T, 



T = T,, 




FIG. 1: Location of the holonomy eigenvalues 
/xi, fi2 shown by (blue) darker small circles, for high 
temperatures (a) and for T — T c (b). Their differ- 
ences v\ — 112 — Mi 5 ^2 = M3 — M2 are also indi- 
cated. The (yellow) lighter circles on the left at the 
phase 7r indicate the case of physical anti periodic 
fermions. 



name 


E 


M 


mass 


M 


± 


± 


V 


M 


+ 




V 


L 






2ttT-v 


L 




± 


2ttT-v 



TABLE I: The electric and magnetic charges and 
the mass (in units of 8tt 2 / g 2 T) for the 4 different 
kinds of SU(2) dyons. 



define a partition function of the dyonic ensem- 
ble which is amenable to practical many-body 
computer simulations, perform the first numer- 
ical simulations and report on the results. 

The partition function of the dyon ensemble 
can be schematically written as 



= J {dXi}e 



- s " det G detF 2 , 



det'F„ 



Vdet'B 

(5) 

where the first factor is the product of the 
differentials of all the collective variables, the 
exponent contains the classical actions of all 
dyons, while the two subsequent determinants 
of matrices G and F zm are related to the 
bosonic and fermionic zero-modes, respectively. 
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The two remaining (primed) determinants cor- 
respond to the nonzero-modes of the one-loop 
bosonic and fermionic operators, describing the 
small linearized perturbations of the classical 
solution. (As usual, the latter determinants are 
divergent and need to be regularized, e.g. in 
the Pauli-Villars scheme, so that the UV bare 
coupling g in the action is combined with the 
UV cut-off into a physical renormalized charge.) 
A quantitative characterization of the dyonic 
ensemble can be achieve by performing Monte 
Carlo sampling of the partition function 

Before we come to such simulations, let 
us outline the physical phenomena to be ad- 
dressed. As the temperature drops from T ~ 
3T C to T c , the value of the holonomy changes 
from zero to its "confining value" v = ttT. 
In the former case the action of the M-dyons 
is small and of the L ones is large, while at 
T < T c their actions become nearly equal. At 
high T "heavy" L dyons can be bound with an- 
tidyons L by the fermion exchanges into certain 
LL clusters [22 , with zero magnetic but dou- 
ble electric Abelian charge. The formation of 
such clusters are reflected in the double-bump 
shape of the Dirac eigenvalue spectrum and a 
gap around zero: in the gas of clusters the chiral 
symmetry remain unbroken. As T decreases, 
the action per dyon decreases and their density 
grow. Also the mass of the fermions decreases. 
Eventually, at some T = T c , the density is high 
enough to bridge the gap in the Dirac eigen- 
value spectrum and produce finite density of 
eigenvalues at zero. This means that the chiral 
symmetry gets spontaneously broken. 

Let us now briefly outline the structure of 
the paper. The next section describes bosonic 
interactions of the dyons, which are sepa- 
rated into approximations made to describe the 
so called moduli space metric associated with 
bosonic zero-modes, and electric screening as- 
sociated with the nonzero bosonic modes. The 
next section describes our approximation for 
the third essential ingredient of the partition 
function, the fermionic determinant. 

The subsequent two sections provide further 
technical details about the partition function 
([5]). Some introduction about the instanton 
dyons and the setting used in our simulations is 



in EI 



reported in section IV The results of the sim- 
ulations are presented and discussed in section 
\V\ while the main conclusions are summarized 



II. BOSONIC DYON INTERACTIONS 
A. Moduli space metric 

In the selfdual (or anti-selfdual) subsectors, 
the dyons are so-called Bogomolny-Prasad- 
Sommerfeld (BPS) protected objects, which 
means that their total action is completely de- 
termined by the total topological charge. Thus, 
at the level of classical solutions of the Yang- 
Mills equations, the total action is independent 
of the dyon collective coordinates. The corre- 
sponding degenerate manyfold of all classical 
solutions with a particular topological charge 
Q, known as the "moduli space", has cer- 
tain metric, which can be calculated from the 
bosonic zero-modes of the classical solution in 
a standard way. The resulting geometry is not 
simple: for monopoles it was derived in classic 
works of Atiyah and Hitchin [25]. Even a con- 
figurations consisting of two monopoles/dyons 
gives raise to a quite nontrivial manifold, with 
a specific singularity known as "bolt" (or ar- 
rowhead). Many-body moduli are extremely 
complicated, and to our knowledge their exact 
metric and geometry has not been explored in 
any detail. 

Obviously, one needs to make some approxi- 
mations, in order to perform practical simula- 
tions. Following Gibbons and Manton [26 and 
Diakonov [21] , the invariant volume element for 
moduli metric can be approximated by the de- 
terminant of a certain matrix G which is de- 
fined in order to incorporate the known limiting 



y/det 



g^detG 



(6) 



The Jacobian determinant of a single L-M 
dyon pair has been calculated by Diakonov, 
Gromov, Petrov and Slizovskiy (DGPS) [24] . 
The nonzero-modes lead to the so called screen- 
ing p henomenon, which we discuss in section 
For the SU(2) gauge group, the matrix G 



II B 



G 



Tr LM Tr LM 

1 A rrr-i I _L 1 



(7) 



Trr 
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where the matrix raws and columns correspond 
to the two dyon types L, M and vlm = \%l — 
xm\- We emphasize that the expression ^ 
turns out to be correct at all distances, not only 
for asymptotically large separations. It was also 
shown that in the limit of trivial holonomy (i.e. 
for fi m — > 0) or at vanishing temperature, the 
resulting measure with ^ reduces to the stan- 
dard 't Hooft single-instanton measure. 

For sake of completeness we report also the 
generalization of Eq. ^ to an arbitrary 
SU(N ) group, 



where 



[G]r 



5 mn (47rv n 
1 



T\x m x m _|_i 



T\x r , 

)- 



Xfn—1 1 



Tlx, 



x m +i\ 



T\x n 



Xm—l | 



(8) 



where the indexes m, n run over the different 
types of dyons. 

The effective potential generated by the mod- 
uli metric can be defined as usual as V e ff = 
— Tlogdet(G). Upon expanding it in powers of 
1/r, one recovers the usual Coulombic potential 
at large distances. 

A classical Coulomb gas is known to be un- 
stable: particles can fall onto each other. For- 
tunately, a pair of non-identical dyons is not af- 
fected by such a problem. Indeed, the Coulomb 
term is complemented by the higher powers of 
1/r of opposite signs. As a result, the total ef- 
fective potential turns out to be only logarith- 
mically divergent at small distances, thus the 
integral which appears in the partition func- 
tion, 



d 3 r det G 



d 3 r 



const. 



(9) 



is well convergent. This expression implies that 
the probability of finding two dyons at the same 
point vanishes. 

Let us now discuss the case in which there 
are dyons of the same kind. Simplifying the 
famous Atiyah-Hitchin many-monopole metrics 
[25] , Gibbons and Manton [26] have shown that 
the weight for K identical dyons reads: 



W =— I dxL-.dxK det G ident , (10) 



ident 



X^ K 2 • _ • 

2-ffc#i=i T\ Xl -x k \ ' 1 - 3 

T\xi-Xj\ 



Note the inversion of the signs here, as com- 
pared to the case of non-identical dyons. We 
emphasize the this expression is exact only for 
large inter-dyon separations, and corrections 
may be important at small distances. 

For a pair of identical dyons, detG = 
Aitv(Atti> — 1/Tr), which remains positive only 
for distances larger than some "core size" r > 
1 / (AttvT) . Note also that the volume element 
detG vanishes linearly at the core, which fol- 
lows from the famous "bolt" geometry. 

Diakonov [21] nicely combined the metric 
tensors for different-kind and same-kind dyons 
into one symmetric-looking metric, describing 
the moduli space of arbitrary number of dyons. 
Let their numbers are K\ dyons of kind 1, K2 
dyons of kind 2,. . . , Kjy dyons of kind N. The 
result is a matrix whose dimension is the total 
number of dyons (K\ + . . . + K N ) x (Ki + . . . + 
K]\f) and reads: 

[^]ra, i,n,j $mn$ij I 47TZ^ m 2 ) — — 

+X] i — ~ 1 + Y1 1 — z~ 1 

^ \%m,i %m-\-l,k\ , \%m,i %m—l,k\ 



%m,k | 



+2 



, n — 1 



\Xrn,i %m,j\ \%m,i %m+l,j\ 
^m,n— 1 



|^m,i ^m— 1, 



(12) 



In this expression x m ^ denotes the coordinate 
of the i th dyon, of kind m. 

Let us now discuss the short-distance be- 
havior. Now, the effective potential V e ff = 
exp(— logG) has multi-particle terms with all 
powers of 1/r. To make sense of it, it is in- 
structive to first calculate it for some examples 
of configurations. For example, we can consider 
a "square" consisting of 2 L and 2 M dyons. 
In is found that, in all cases, their combined 
effect can be described by a weakening of the 
Coulomb divergence. 

When the distances between dyons become 
too small, the Coulombic term becomes unrea- 
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sonably large, and it needs to be regulated: af- 
ter all, a dyon is not a point charge and has a 
certain size. In the present work, we adopt a 
regulation based on the substitution 

o o 

(13) 



r vV 2 + a 2 ' 

where we have chosen in the simulation the 
value of the cutoff parameter a = n/T, leading 
to vanishing measure for two identical dyons, 
in the limit r — >• (to be referred to as "soft 
core" ). This is in contrast with imposing that 
the measure should vanish at some finite r, as 
in the famous Attyah- Hit chin metric (to be re- 
ferred to as "hard core"). 

The configuration space of a multi-dyon 
ensemble can be sampled with a stan- 
dard Metropolis-based Monte Carlo algorithm. 
Within such an approach, we shall not only en- 
force the global positivity of the metric deter- 
minant, but also we will require each eigenvalue 
of G to be positive. This condition is becom- 
ing restrictive at high dyon densities, i.e. when 
fraction of the total volume filled by dyons is 
VT 3 ~ O(l), eventually making the Monte 
Carlo sampling of a very dense dyon ensemble 
rather inefficient. (This kind of computational 
difficulty is well known, e.g. in simulations of 
classical liquids with hard repulsive cores.) 



B. Nonzero- modes and electric screening 

The DGPS determinant of the non-zero- 
modes [24] for a single LM dyon pair contains 
the so-called "confining" term, proportional to 
the inter-dyon separation and to the "screening 
mass" 



logdet'i? = tlm x 



2ttM 2 



which displays the so-called electric Debye mass 
M 2 = (2/3)<? 2 T 2 , in the SU(2) gluondynam- 
ics. Upon generalizing to any N c and number 
of fundamental fermions Nf this expression be- 
comes [27] . 



M 2 D 



g 2 T\N c /3- 



■N f /6). 



(15) 



We now briefly elucidate the origin of the 
linear potential which appears in Eq. (14). 



Let us start with the quartic term [A^A^] 2 of 
the Yang-Mills Lagrangian, and imagine two of 
those fields are classical ones, say from some 
dyon solution, and two others belong to thermal 
perturbative gluons of the heat-bath. While the 
latter produce the Debye mass, the former gen- 
erate the linear dependence on r^, from the 
volume average of the squared Coulomb large- 
distance potentials 



V 12 



d 6 x 



1 



1 



4tt r LM 



(16) 

Note that we consider here neutral L — M pair 
(i.e. the splitting of one instanton), as a result 
of which their Coulomb potential cancels out at 
large distances from the pair and the volume in- 
tegral is convergent. Clearly, non- neutral con- 
figurations cannot be treated this way. 



We further note that the form ( 14) can be ob- 



tained directly by the instanton screening term 
calculated by Pisarski and Yaffe [28] by recall- 
ing that the instanton size p and the L — M 
separation are related by the expression 



irp T = tml- 



(17) 



which relates the "4-d dipole" of the instanton 
field to the "3-d dipole" of the dyon LM pair 
made of opposite charges. 

Let us now work out the corresponding gen- 
eral formula for screening potential which holds 
in the many-body case. The sum over all dyonic 
contributions to A^ can be written as 



I 



d 3 a 



(18) 



where now the sum runs over all dyons with 



(14) Qi = =1=1 is the charge and = \x — Zj 



One can write (A4) 2 as a double sum, in 
which we separate the diagonal terms i = j 
from all non-diagonal terms. Unless total neu- 
trality is ensured 



^2Qi = o 



(19) 



this integral is divergent at large distances. If 
there is an overall neutrality, one can regulate 
the sum term by term, by subtracting the cor- 
responding (r-independent) divergency. Let us 
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say discuss a non-diagonal term i ^ j which we 
can then be rewritten as 



2 



— >> 



1 



1 



'i'j 




reducing it to the 2-body case above ([16]). 
Therefore the total answer is a simple gener- 
alization of the linear potential 



(21) 



i>j 



with Qi = ±1 being the electric charges. In 
contrast to linear potential induced by the con- 
fining flux tubes, in this case there is a sum 
over all pairs of dyons. The sum therefore has 
N(N — 1)/2 terms for N dyons, of various signs, 
which together enforce local neutrality of the 
dyon ensemble, making as small as possi- 
ble. It does not diverge in infinite volume be- 
cause distant dyons may be combined with their 
counterparts of the opposite charge, producing 
dipole fields. On a manifold with the topology 
of a sphere extra care is needed, as discussed in 
Appendix [B] 



III. THE FERMIONIC DETERMINANT 

In evaluating the fermionic determinant, the 
Dirac operator is approximated by retaining 
only the contribution evaluated on the sub- 
space of fermionic zero-modes of the individual 
pseudo-particles |0q): 

Det(i7 /i D /i + im) ~ Det(f + im), (22) 
where 



(23) 



This scheme was well tested in the frame- 
work of the instanton liquid model, where 
it corresponds to summing up all loop dia- 
grams created by 't Hooft effective Lagrangian. 
One of the most important lessons from that 
was existence of rather narrow "zero-mode 
zone" (ZMZ) of the Dirac eigenvalues |A| < 
a ~ 20 MeV, which drive all the effects related 
to chiral symmetry breaking. The narrowness 
followed from relative diluteness of the instan- 
ton ensemble. 



The dyonic formulation also generates a nar- 
row ZMZ: in fact, at nonzero holonomy, all light 
quarks become effectively massive, so the so 
called "hopping amplitudes" - i.e. the non- 
diagonal matrix elements of the Dirac opera- 
tor between different dyons - are exponentially 
suppressed with distance, 



T 



exp(— M r) 



(24) 



where r is the inter-dyon distance. M is the so- 
called holonomy mass M, should not be con- 
fused with the ordinary (i.e. current) quark 
mass m. In particular, the former does not 
break chiral symmetry, since it arises from the 
interaction of the quarks with the A4 compo- 
nent of the gluon field. An estimate of the mag- 
nitude of M for various fermions followed from 
the zero-mode solution can be found e.g. in Ref 
[22] . Antiperiodic fermions in SU(2) have zero 
modes on L dyons, with the (holonomy-induced 
chirally symmetric) mass given by 



M = ttT-v/2 



(25) 



proportional to the distance from the fermion 
boundary condition angle ( the (yellow) lighter 
circles on the left at the phase ir in Fig|T]) 
to the (blue) dark circles indicated the holon- 
omy values. If the fermionic boundary phase is 
changed, the yellow circle moves and then one 
should select the smallest distance between the 
two holonomy values. 

The most important qualitative trend is the 
following: as the temperature decreases from 
high T to T « T c , the holonomy changes as 
indicated in Fig (I] As a result, the effective 
fermion mass M decreases and the "hopping 
amplitude" p4| ) increases. Another qualitative 
tendency is of course the reduction of the effec- 
tive coupling in the dyonic action, which addi- 
tionaly increases the dyon density and reduces 
the inter-dyon distance r. Both effects eventu- 
ally conspire towards breaking the chiral sym- 
metry. 

Some of the (block) entries in Eq. (23) 
are identically vanishing. In particular, the 
anti-periodic conditions on fermions imply that 
zero-mode exist only for L and L type dyons, 
see e.g. [22 for explicit solution. Hence, we 
can ignore the overlap zones involving M— and 
M— type dyons. In addition, the Dirac oper- 
ator mixes only modes with opposite chirality, 
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hence also Tll 
operator reduces to: 



(i j^D^ + im) 



0. Hence, the Dirac 



(26) 




where the the ij element in LL block is given 
by the (approximate) formula [22] containing 
the "holonomy mass" 



-Mri 



LL 



V 1 + Mr ij 



(27) 



Notice that the normalization constant c is ir- 
relevant for Monte Carlo ensemble simulations: 
it may only be needed if some physical signif- 
icance will be ascribed to the quark mass m 
rather than pushing it to zero. For Nf massless 
fermion flavors the determinant simply reads 



det it) = | det T lL 



27V f 



(28) 



A. The case of a single dyon molecule 

Let us first consider the case in which there is 
only one dyonic "molecule" made by L-, Z-, M- 
and M- type dyons. The fermionic determinant 
can be viewed as the second-order one-loop dia- 
gram in 't Hooft effective Lagrangian, with two 
vertices with 2Nf fermionic propagators in be- 
tween. From Eq. ([28]) and p7|) we get 



det iD 



\2N f 



l LL\ 



-2N f Mr LL -\og(l+Mr LL ) 
(29) 

Neglecting the logarithmic dependence, the 
formula for the determinant for 1 molecule 
(L, M, Z, M) reads as en effective potential 



V = - log det iD 2N f M r LL , 



(30) 



where again we have dropped an irrelevant nor- 
malization factor. Thus, fermion exchange also 
creates a linear confining potential, but now be- 
tween dyons and antidyons. Also, note that 
there is no minus sign and that an additional 
parameter - the number of flavors Nf - ap- 
pears as a factor. Therefore, one expects that 
fermion exchange to generate the tightly bound 
LL clusters first described in Ref 22 . 



B. The case of two dyon molecules 

This is the simplest example demonstrating 
how dyonic clusters undergo mutual repulsion. 
To see this we study the two dyon molecule 
case, restricting to the L-dyon zone and label- 
ing the pseudo particles with Li, L2, L\ and L2 
we find: 



i D 



Hence 



( * 




L 2 


u 


U 


Li 








T L 1 L 1 


T L 1 L 2 


L 2 








T L 2 L 1 


Tl 2 L 2 


U 


T f - 


T t _ 

L2L1 








\ L2 


T f - 

L\L<2 


T f - 

L2L2 









(31) 



det iD = -|detT i£ | 2JV ' 



-\T Ll - Ll T:> 



L 2 L 2 
2N f 



(32) 



from which one immediately sees that the 
configuration weight vanishes when the two 
molecules overlap. In the sext sections we shall 
we shall perform numerical simulations which 
will allow us to investigate the consequences of 
such a "Fermi repulsion". 



IV. THE SETTING 
A. The setting and parameters 

Let us start by enumerating parameters of 
our physics problem, which we will group into 
global, internal and external ones. 

• The global parameters are the number of 
colors N c of the theory and the number 
of fundamental fermions Nf of degener- 
ate current mass m. In this work N c = 2 
and Nf = 0, 1, 2, 4. One may in addition 
consider other types of fermions, e.g. ad- 
joint ones, but we leave such an extension 
to future work. 

• The internal parameters are those on 
which the partition function depends, 
such as the heat-bath temperature T, and 
the flavor chemical potentials \±i. Al- 
though the N c = 2 theory allows for fi- 
nite density Monte Carlo simulations, for 
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now we will only discuss the zero density 
case, Hi = 0, using T as the single internal 
parameter. 

• The external parameters are those which 
may be imposed. For example, one 
may consider modifications of the QCD 
vacuum in the nonzero (QED) magnetic 
field: it is now in fact a quite active re- 
search field. Another important example 
of an external parameter is the CP-odd 
angle which however have a sign prob- 
lem and need reweighing or analytic con- 
tinuation. None of it will be discussed in 
this work. 

The dyonic model is defined by the partition 
function defined above and can be studied for 
any values of such parameters. After this is 
done, the results can be mapped/compared to 
lattice data using appropriate values of these 
parameters as inferred from lattice simulations 
at the corresponding temperature. 

One such input is the value of the holon- 
omy parameter fi at temperature T. There are 
rather accurate lattice data defining it, see Ap- 
pendix [Cj 

Other simulation inputs are the densities rii 
corresponding to dyons of various kind i. We 
will also use their proper volumes defined by 
their inverse densities, i.e. Vi = l/ft*. At the 
time of this writing the lattice studies of the 
dyons are too few to allow for a systematic map- 
ping, of the dyon densities at different temper- 
ature. While we do not really use such data, for 
the orientation of the reader we put brief review 
of the information we find in the Appendix [D] 

We will only mention here that since the 
masses of L and M dyons are different, their 
densities on the lattice are different as well. 
Note that such disparity between the densities 
of the L and M dyons produces a nonzero den- 
sity of the Abelian electric (but not magnetic) 
charge. This is of course not a problem in gen- 
eral, as it is compensated by an asymmetry of 
the densities of the electrically charged gluons. 

It would however create a technical problem 
for our simulation setting, as we focus entirely 
on the topological (dyonic) sector. Further- 
more, one has to use a compact many fold - a 
torus or a sphere - in which all the charges must 
be compensated and field lines accounted for: 
hence only totally neutral systems are allowed. 



One possible solution to this problem is simply 
adding to each charge the homogeneously dis- 
tributed compensating charge: see Appendix [B] 
for details. Another practical solution we adopt 
in this work is to ignore the extra M dyons: so 
we simulate an equal number of the M, L dyons. 

Finally, about our units. Standard nonper- 
turbative normalization of the gauge field elim- 
inate the coupling from the Yang-Mills equa- 
tions. It appears in classical action S = 
0(1/ g 2 ), but since its total value depends on 
the total number of dyons which is not changed 
in our configurations, it is absorbed in the over- 
all input density. The one-loop bosonic and 
fermionic determinants we actually simulate do 
not contain the coupling. Thus, there is no 
running and scale dependence in the model at 
all, and all quantities are essentially dimension- 
less. Or, as we do, all quantities can always be 
made dimensionless by multiplying by appro- 
priate powers of temperature T. 

In particular, the simulations are done for 
fixed number of the dyons on a S 3 sphere, with 
its radius characterized by the (dimensionless) 
radius R(S 3 )T, or by the proper volume per 
dyon in T-units, VT 3 , the inverse of the di- 
luteness (we remind that the total volume is 
V(S 3 ) = 2tt 2 R(S 3 ) 3 . The values for which the 
simulations has been done are listed in the Ta- 
ble |llj Note that the density or volume per dyon 
changes by about two orders of magnitude. The 
largest values of the proper volume (per dyon) 
correspond in physical ensemble to weaker cou- 
pling, larger dyon masses and higher temper- 
atures. Going down in this volume one goes 
into the denser systems, stronger coupling and 
the lower T. Respectively we will monitor in 
simulations how all observables depend on the 
proper volume, focusing on the correlations be- 
tween the dyons and the dirac eigenvalues spec- 
trum. 

The only dimensional parameter which en- 
ters the fermionic determinant is the current 
quark mass m, which in our units corresponds 
to the ratio m/T. Keeping its constant means 
a temperature-dependent mass. This does not 
represent a problem as long as one is interested 
in the massless limit, m — >• 0, in which the 
classical action becomes scale invariant. Yet in 
reality, following a prescription which is com- 
monly adopted by lattice practitioners, the ac- 
tual value of this mass used in our simulations 
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FIG. 2: The correlation function for LM,LL and LL dyons versus distance, normalized to the volume 
available. From top to bottom we show Nf = 1, 2,4, respectively. Left/right columns are for the volumes 
per dyon VT 3 0.31, 1.04. 



is not zero and is not directly related to physical 
light quark masses: its role is to prevent influ- 
ence of the finite volume effects. In particular, 
we shall will see below that, in order to cor- 



rectly monitor the chiral symmetry breaking, 
we are bound to consider masses m/T > 0.2, 
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dyons 


R(S 3 )T 


VT 3 /dyon 


64 


4.5 


28. 


64 


3.0 


8.3 


64 


2.5 


4.8 


64 


2.2 


3.28 


64 


1.5 


1.04 


64 


1.2 


0.53 


64 


1. 


0.31 



TABLE II: The list of the radii and volumes/dyon 
used in the simulations 

V. SIMULATION RESULTS 

One standard observable in the simulation is 
the acceptance probability of Metropolis steps, 
which should be tuned to being comparable to 
rejection probability to run it efficiently. We es- 
timate the number of Metropolis steps needed 
to achieve thermalization by standard auto- 
correlation analysis based on the dyons' posi- 
tion. Our typical runs include a total of 256 
independent configurations, obtained from 32 
independent Markov chains, each consisting of 
64000 Metropolis steps. All simulations were 
performed at the Interdisciplinary Laboratory 
for Computational Science (LISC), at Trento. 

A. The spatial correlations between dyons 

One of the benefits of going from 4-d torus 
to the S 1 x S 3 geometry is that in the latter 
case it is much easier to define the interparticle 
distance 

rij = R aij = R arcos (n^nj) (33) 

where cos(Sij) is the angle between two points 
defined via the scalar product of their unit po- 
sition vectors on S 3 and R is its radius. 

We collected histograms of the various two- 
particle correlation functions as a function of 
the (dimensionless) distance RT on the sphere 
are normalized by dividing out the volume ele- 
ment dV/dSij ~ sin 2 (5ij) corresponding to ran- 
dom occupation on the S 3 sphere. A sample of 
the results for different values of Nf and differ- 
ent dyonic densities is shown in Fig. [2] One can 
immediately identify a number of feature which 
are common to all such correlation functions: 




K/T 



FIG. 3: The spectrum of the Dirac eigenvalues in 
the quenched ensemble, comparing two ensembles 
of 64 and 128 dyons at the same density, with zero 
current mass. The shaded area denotes the region 
of the spectrum where volume artifacts become sig- 
nificant. 



• repulsive correlation between identical 
dyons, 



• attractive interaction between the dyon 
and antidyon; 



• the radius of the correlation decreases as 
Nf grows 



B. The Dirac eigenvalues and the final 
size effects 

Our present computer resources allowed us 
to perform simulations with 64 dyons. In or- 
der to estimate the magnitude of finite- volume 
effects, we performed also a few quenched runs 
with ensembles consisting of 128 dyons. Two 
distributions of the Dirac eigenvalues for 64 and 
128 dyons performed at the same density con- 
ditions are compared in Figj3J One can see that 
the finite volume effects are significant only for 
the smallest eigenvalues. More specifically, one 
finds distortions of the spectrum for the first 
bins, say for the shaded area 

Y < 0-25 (34) 
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Notice also that that larger systems have 
smaller finite size effects, as expected. 

A rather crude account of the finite-volume 
effects can be thus made by simply disregarding 
the data points in the strip identified by the 
inequality (34), yet we will need more refined 



fits when we introduce a nonzero quark mass in 
the next section. 



C. Quark condensate dependence on the 
current quark mass 

Let us start reminding that there are two 
types of quark masses: "current mass" gen- 
erated by the interaction with the Standard 
Model scalar (Higgs) VEV with the left-right 
structure 



L L r = m (q L qR + qRqL) 



(35) 



and thus violates the chiral symmetry explic- 
itly. The other "mass" M is generated by the 
holonomy. As it is due to the VEV of the vec- 
tor field component A±, its structure is chirally 
diagonal (symmetric) 



La a = M(q L q L + qRqRj 



(36) 



We recall that M appears in the parametriza- 
tion of matrix elements of the T-matrix, which 
defines the Dirac operator in the background 
holonomy potential. (For clarity, we do not use 
the notation of "constituent quark mass" in- 
duced by sponateous chiral symmetry breaking 
in this paper.) 

The current quark mass m is important for 
our study because it "seeds" the chiral symme- 
try breaking. As it is well known to lattice prac- 
titioners, in order to study this spontaneously 
broken symmetry one formally should take the 
double limit, the infinite volume V —> oo and 
massless quarks m — » 0, in this order. In prac- 
tice, all simulations are performed on a finite 
volume, and the absolute values of the smallest 
Dirac eigenvalues are bounded from below at 
the scale A ~ 1/V. The current mass should be 
selected somewhat above this scale, in order to 
cause the chiral symmetry breaking. The rea- 
son is obvious from the expression for the quark 
condensate, as the trace of the massive quark 



propagator 



m = E 



m 



Xf + m 2 



rn I dX ,/ (A) 2 (37) 



Tuning the quark mass values used in the sim- 
ulation is important: it should be as small as 
possible, yet large enough to avoid large finite- 
volume effects. 

An example of how this tuning works in prac- 
tice is demonstrated in Figj4j In the left panel 
we show the Dirac eigenvalue distributions for 
6 values of the current quark mass. As one 
can see, the spectrum itself changes, from a dip 
near A = at zero mass, to a sharp peak at 
larger masses (as for quenched Nf — case). 
The Casher-Banks "density of the eigenvalues 
at the origin" thus changes dramatically, from 
zero to very large values: hence, reading off p(0) 
is clearly not a good method for estimating the 
quark condensate estimate. By contrast, the 
points shown in the right panel of Figj4] show 
the quark condensate estimated from the loop 
expression (37). We observe that the quark 



condensate remains approximately mass inde- 
pendent for larger masses m > 0.2T, while de- 
creasing for smaller masses, eventually going to 
zero at m = 0. 

Let us first comment on the continuous lines 
in FigQb), passing through our 6 simulation 
points. Those correspond to the loop expres- 
sion ( |37| ) in which one uses the eigenvalues as 
determined from the simulation with the corre- 
sponding dynamical mass, but in the loop for- 
mula itself the mass is substituted by a separate 
parameter, called the "valence quark mass" m v 
in literature. Such "valence approximation" 
has been used before, sometimes for quenched 
ensembles, to approximate the mass depen- 
dence outside of the region where simulations 
were done. However, as seen from the plot, the 
trend indicated by those curves is quite mis- 
leading, being completely different from what 
is obtained from the dynamical simulations - 
the points themselves. It shows that the main 
mass dependence is stored in the eigenvalue dis- 
tribution, reflecting the dynamical quark mass 
of the vacuum, of the "sea" quarks. 

Let us now remind the Smilga- Stern theo- 
rem [31], which is based on chiral perturbation 
theory and predicts the slope of the eigenvalue 
spectrum near the origin, as a function of the 
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FIG. 4: Left panel: Dirac spectrum as a function of the current quark mass m. Ri ght panel: quark 
condensate as a function of the current quark mass obtained from the loop-estimate p7| . The circles 
denotes the value of the mass used in the generation of the dyonic configurations. The simulations used 
for the results shown in both panels were performed on a S3 of radius 1.5T and with Nf — 2. 



chiral condensate (qq) , the pion decay constant 
and the number of flavors: 



dN 



(A) 



1 



|A| 



dX 

4 ) (QQ) 2 



(qq) (38) 



N 



f 



327r 2 Fi 



0(A 2 ) 



Lattice and instanton liquid simulations had in- 
deed confirmed that the slope of the distribu- 
tion changes sign at Nf = 2, from a negative 
constant at Nf < 2 to positive at Nf > 2. 

Let us also point out that the spectral den- 
sity in the expanded form near A = can be 
put into the integral in (37), generating the 



condensate expanded in powers of m. The ze- 
roth order is the Casher-Banks term. Return- 
ing to Figj4|b) for Nf = 2, we see that no liner 
term in m is expected, and higher order terms 
can be neglected: so the correct infinite volume 
V — >• 00 limit should be just a constant. Our 
guess of its value is indicated by a horizontal 
dashed line. Comparing it to the data we see 
that for m/T < 0.2 finite solum effects are large 
and those should not be used. Since we cannot 
afford to do simulations with variable masses 
in all cases studied, we will below only use the 
dynamical mass m/T = 0.25, at which we ex- 
pect the condensate accuracy determination to 
be not worse than say 20%. 

Our computed eigenvalue spectra, for differ- 
ent volumes and number of flavors, are shown 
m Figj5] While some finite- volume effect is still 



seen at the smallest eigenvalues, it is much less 
prominent than in zero mass simulations dis- 
cussed above: they are only seen for X/T < 0.1. 

The Nf = 4 spectra are clearly "gapped", 
without any small eigenvalues, and the gap 
magnitude grows as the dyon density falls. We 
thus conclude that Nf = 4 case needs rather 
large dyon density to break the chiral symme- 
try, perhaps close to the highest one we used. 

Other Nf simulation show no gap and thus 
correspond to chirally broken phase, except for 
very low dyon density. At different dyon den- 
sity and number of flavor, we performed a linear 
fit p(X) = Q + CX. The fit was restricted to a 
range with 0.1 < X/T until the maximum A 
for which the reduced x 2 remains smaller than 
1. This way, we identified constants Q and 
C with the quark condensate and the Smilga- 
Stern constant, respectively. The result of the 
fit is shown in FigjBJ The nonzero conden- 
sate observed for the Nf = 1 and Nf = 2 en- 
sembles at high dyon density is rather density- 
independent. 

More systematic account for the finite- 
volume effects and more accurate determina- 
tion of the condensate can be done by compar- 
ing our data with the expectations for the so 
called "mesoscopic regime", in which the vol- 
ume is not macroscopically large. Quantitative 
predictions of the shape are known from chi- 
ral random matrix models [33], which are well 
confirmed on the lattice. We plan to do so else- 
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FIG. 5: The Dirac eigenvalue distributions for three values of the Nf = 1,2,4. The plots, ordered from 
top-left to right-bottom are performed on different radius of the 3-sphere (R=1.0 T,R=1.5 T ,R=2.2 T, 
R=3.0 T ,R=4.5 T), from high to low dyon density. The current quark mass was set to m/T — 0.25. 



where. 



D. Free energy, dyons' interaction, the 
back-reaction and confinement 

The well known perturbative holonomy po- 
tential, which has a minimum at zero holon- 
omy, has been argued by Polyakov's original 
work to get cancelled by the nonperturbative 



effects, resulting in a vanishing Polyakov line 
and thus confinement. Lattice study of the ef- 
fective holonomy potential, in particularly re- 
cent Ref.[34 , have indeed found such behavior. 
Diakonov [2Tj further argued that the nonper- 
turbative contribution comes from the back re- 
action of the dyons. He also argued that (at 
least the Coulomb-like moduli part of) the dyon 
interaction is small and can be approximately 
ignored. In this section we make the first step 
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FIG. 6: The quark condensate versus the dyon 
density. The current quark mass is set to m = 
0.25 T. 



FIG. 7: The dependence of the average interaction 
free energy on the dimensionless adiabatic param- 
eter A. 



toward understanding of whether those ideas 
are correct. 

We calculate the change in the free energy of 
the system, between the interacting and nonin- 
teracting dyon gas using standard a thermody- 
namic integration method, based on the adia- 
batic switching of the interaction between the 
particles. One splits the action into indepen- 
dent particles and their interaction, introducing 
the adiabatic parameter A 



S — Sq H~ AS^ 



int 



(39) 



and do simulations with A changing from to 1. 
The resulting free-energy is recovered as follows 



F = Fn 



f 

Jo 



dX < Si n t > 



(40) 



The resulting dependence of the action as a 
function of the adiabatic parameter is shown 
in Fig(7| Of course, we include all three ingre- 
dients of the interaction, the moduli, screening 
and the fermions, and measure them separately 
as well: but for brevity we will not discuss those 
details here. 

As one simulates an ensemble with a fixed 
number of dyons, one can disregard any con- 
stant factors in the partition function, as the 
relative weights of the ensemble configurations 
depends only on terms which are functions of 
the dyonic collective coordinates. This explains 
why one has positive free-energy of the nonin- 
teracting ensemble at A = 0, which is not really 
physical and depends of what factors we do or 



do not include in the free gas expression. The 
free-energy change AF = —22 is physical, and 
a shift to negative is typical for liquids. The 
value itself is for the entire system of 64 dyons: 
it thus corresponds to a free energy change per 
particle of AF/N = —0.34 which is due to the 
attraction in the system. 

To put it in perspective, one should know the 
action per dyon. We will discuss available lat- 
tice data in the Appendix: the effective action 
fitted to those are about S ~ 3 per dyon. The 
result of this section show that the average in- 
teraction part of it is reasonably small, ~ 1/10. 
On the other hand, the partition function is en- 
hanced by about exp(.34) = 1.4 per dyon, not 
a negligible enhancement. 

As for the back reaction to the total holon- 
omy potential, we think that it is premature to 
discuss it now, before lattice simulations pro- 
vide more define estimates on the dyon density. 



VI. SUMMARY AND FURTHER 
DEVELOPMENTS 



This paper is the first study of the topol- 
ogy, in the finite-temperature QCD near T c , 
which incorporates the nonzero holonomy and 
the (anti)self dual dyons as the ingredients. We 
(i) formulate the partition function for (anti) 
self-dual dyons incorporating one-loop screen- 
ing, moduli space measure and sermonic zero 
modes; and (ii) perform the first direct sim- 
ulation of the resulting statistical ensembles. 
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We obtain certain set of results, showing how 
the dyon system depends on the dyon density 
(changed by about 2 orders of magnitude), on 
the number of light quark flavors Nf, from 
to 4, and on the (current) quark mass m. We 
have in particularly studied the Dirac eigen- 
value spectra and chiral properties: the gaps 
in the symmetric phase and the quark conden- 
sate in the broken phase. More generally, we 
may say that in spite of the long-range na- 
ture of Coulomb forces involved and linearly 
rising screening, the simulations of the dyonic 
system turned out to be not as challenging as 
we thought when we started this work. While 
the correlations between constituents are lo- 
cally strong, overall liquid-like behavior is ob- 
served. 

At the moment we are not calculating the 
fermionic propagators and correlation functions 
of operators with various quantum numbers, 
but it can be done rather directly, in a way 
similar to that in the instanton liquid. 

Mapping those results to fmite-T QCD can 
be done provided the information about the 
dyons, such as their density we reviews briefly 
in the Appendix, are extracted from the lattice. 
One should also directly relate fermionic quasi 
zero modes with the underlying topology, which 
we believe will display the dyonic dominance 
in the chiral symmetry breaking phenomenon, 
in the near-T c region. One may also expect 
progress in respect with the back reaction of 
the dyons on the holonomy potential and the 
origin of confinement. 

Another direction is to extend the model to- 
ward theories in which one has high-T con- 
finement under the analytical control, see e.g. 
[35] , by including adjoint fermions with modi- 
fied (e.g. periodic or with an arbitrary phase) 
boundary conditions on the torus. 

Extending our results to lower temperatures 
- or much smaller holonomy values - remains 
a numerically challenging task. Attractive 
correlations of the LM pairs should mostly 
recover the original instantons, yet preserving 
small Polyakov line value and confinement. 
At the same time, fermions induce the LL 
attractive interaction, producing topologically 
neutral clusters whose role in the QCD vacuum 
has not been elucidated or studied. 
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Appendix A: Distances on the S3 sphere 

In our program, all distances are defined on 
a 3-dimensional sphere embedded in 4 dimen- 
sions. Let z l n be a unitary vector in 4 dimen- 
sion, which will be associated to the position 
of the pseudo-particle of type i in the n— th 
molecule. Monte Carlo trial moves are made 
by performing random rotations of given small 
angle (boldness) on the sphere. 

The angle between the 4-D vector z % n identify- 
ing the pseudo-particle of type i relative to the 
molecule n and z 3 m that points to the pseudo- 
particle of type j relative to the molecule m on 
the unitary sphere is given by: 

ainjm arccos [4 • z 3 m ] (Al) 

(Notice that this formula is ill-defined for i = 
j D n, m. On the other hand, the distance of a 
particle from itself is never needed in the cal- 
culation.) Hence the distance between these 
pseudo particle is simply: 

Tin,jm — & ^in y jmi (^^) 

where a is a length scale which can be identified 
with the box size. 

For illustration purposes we mapped the S' 3 
to a flat R 3 space using the stereographic pro- 
jection. The three dimensional coordinates of 
our dyons are defined on a sphere, hence we use 
4-dim vectors with the constraint + 7?! +n§ + 
n\ = R 2 . 

One may transform them into the 3 polar an- 
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gle angles like this 



77,1 = R COs(lp) 

ri2 — R sin( , 0) cos(#) 

ns = R sm(ip) sin(0) cos(0) 

= R sin( , 0) sin(0) sin(0) 



(A3) 



The inverse formulas are: 



= cotg 



2 cotg 



ip = cotg 1 



n 2 



n 3 



714 



ni 



(A4) 



In a stereographic projection to the flat 3d- 
space in polar coordinates r, (/> the 9 and <p 
angle keep the same meaning as above, while 
the distance from the origin r is defined by the 
angle ip according to 



r = 2 i?tan(7r/2 - ip/2) 



(A5) 



(R=l is the sphere radius and the argument of 
tan is the angle conjugated to psi in the equi- 
lateral triangle of the projection). 

Note that ip = tt maps the south pole to 
the origin r = 0, and the equator ip = ir/2 
is mapped to r = 2R. Hence the cartesian co- 
ordinates on the projected flat space are: 

x 1 = 2 R tan(7r/2-^/2) sin(0) cos(» 
x 2 = 2 R tan(7r/2 - tp/2) sin((9) sin(0) 
x 3 = x x = 2 R tan(7r/2 - ip/2) cos((9)(A6) 



For a point charge at the north pole (ip = 0) the 
solution is obviously independent on <p and 
we need to solve only the ip -dependent part 



/"(^) + 2/ / (^)/tan(^)=4^ 



(B3) 



For zero r.h.s. one get solution / = 
(l/4ntan('0) which looks like a positive charge 
at one pole and the negative charge at the 
other. Indeed, on a compact manyfold all field 
lines have to go somewhere, and symmetry re- 
quire it to be the opposite pole. 

One remedy can be introducing the homoge- 
neous compensating charge. This solution is 
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1 



tan(^) 



(B4) 



It is compared with the naive Coulomb-like l/ip 
and solution given above in Fig. One can see 
from it that its behavior is very reasonable, and 
it vanishes at the opposite pole. 

Screening effect contains the averaging of the 
dipole potential over the 3-d space, which pro- 
duces the linear rising effective potential (as ex- 
plained in the text). Keeping one charge at the 
north pole and spreading another at ip = a 2d 
sphere one gets another solution of the Poisson 
eqn. It is 

A (& <a) = 1 W" - i W" (B5) 

v ' 4 tan(^)7r 4 tan(a)7r v ' 

with A(ip > a) = 0. Integrating the square of 
this field over the 3d sphere one gets the follow- 
ing expression 

f a ii a2 • / ,\2 a-sin(a) cos(a) 
/ #A 2 sin = 2 V / v 

Both agree reasonably well till about the equa- 
tor of the 3d sphere. 



(B6) 



Appendix B: The Coulomb fields on a S 3 

The usual coordinates ip,6,(p have the follow- 
ing line element (metric) 

dl 2 = dtp 2 + sin(iP) 2 d6 2 + sin(iP) 2 sm(0) 2 d(P 2 

(Bl) 

Standard Poisson equation is 



1 

71 



dpy/gg^dvf = 4n P (B2) 



Appendix C: SU(2) holonomy on the lattice 

One important potential input is the holon- 
omy value v(T) (i.e. the VEV of the Polyakov 
line's VEV (P)). While we have not yet used 
it, we indicate that its temperature dependence 
for the SU(2) gauge group is known. To our 
knowledge the most accurately determined val- 
ues are from Ref. [32]. The resulting fit of the 
Polyakov line VEV 

(P(T)) = 1.22 t - 326 (1 - 0.246 t - 53 ), (CI) 
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where t is the near-critical parameter 

t = (T-T c )/T c (C2) 

The temperature dependence of the holonomy 
is then arccos(P) = v(T)/2T. 

Note that these data and the fit are consis- 
tent with the indices of the Ising universality 
class to which the SU(2) pure gauge theory be- 
longs. 

Appendix D: The (anti)self-dual dyons on 
the lattice 

While we only consider a case with 2 colors, 
we remind the reader that one of the major ben- 
efits of the dyon approach is the smooth behav- 
ior of the action in the large- N c limit. Indeed, 
the dyon partition function is roughly the 1/N C 
root of the instanton partition function, so the 
dyon action in the large- N c limit is expressed 
via 't Hooft coupling A = g 2 N c as 



While at very strong coupling A — > oo it is 
small, at the "physical" value of the parameter 
A = 20 — 30 usually used for strongly coupled 
QGP applications it is still in the range of 3-4 fir- 
units. Hence, the weight factor exp(— S) is still 
a rather small number, providing a substantial 
suppression of the dyon density. It is important 
to stress that, although in principle it may be 
argued that the value of the action So ~ 3 — 4 
is not large enough to justify the usage of the 
semiclassical approach, the corresponding dyon 
ensemble remains dilute even close to T c , hence 
providing an a posteriori justification of the ap- 
proach. 

While we generally postpone mapping of the 
results of our dyonic model to the QCD ther- 
modynamics directly, in this section we still dis- 
cuss some available results on the dyons, from 
Refs. [29] and [30]. One of the issues here is 
whether they are or are not compatible with 
the with semiclassical expressions, in spite of 
possible low accuracy of those. 

These works both investigate the SU(2) glu- 
odynamics on the 20 2 x 6 and 20 2 x 4 lat- 
tices respectively, adopting the same tree-level 
improved Symanzik action at a bare coupling 



/3 C = 2N c /g 2 = 3.248. In the former lattice, 
this choice corresponds to the critical temper- 
ature, T = T c , while in the latter lattice to a 
temperature T = (3/2)T c . 

Note that lattice practitioners often intro- 
duce ad-hoc units of energy and length which 
are defined in such a way to mimic the dynam- 
ics of the physical world. For example, In order 
to fix an energy scale, they measure the dimen- 
sionless constant. 

\/cr(T = 0)a = 0.236, T c /^cr(T = 0) = 0.71 

(D2) 

and artificially set the vacuum string tension 
a(T = 0) value to (420 MeV) 2 , i.e. to its 
physical value in the 'real- world' QCD. Fur- 
thermore, one keeps also the 'real-world' rela- 
tion between the length and the energy units 
fm = 1/(197 MeV). Under this set of conven- 
tions, one obtains the following "absolute" val- 
ues for the lattice spacing and the critical tem- 
perature 

a = 0.236 x (197/420) ~ O.llfm (D3) 
T c = 0.71 x 420 « 300MeV. (D4) 

When comparing to lattice results from those 
works - such as the density of various dyons - 
we may adopt such absolute units. However, 
in general, we will refrain from doing so be- 
cause this may lead to confusion when study- 
ing the dependence on the number of flavors 
Nf. In addition, holding fixed the value of the 
string tension constant is only practical in a lat- 
tice approach, in which this observable can be 
straightforwardly calculated. In view of such 
considerations, in the present analysis we avoid 
using "physical units" and consider only dimen- 
sionless ratios and couplings. 

In Ref.[29 , Ilgenfritz and co-workers report a 
list of the absolute densities of static (M-type) 
dyons and calorons (LM dyon pairs) which 
were found in their lattice calculation. As M 
dyons are much lighter than L, the calorons 
density is dominated by the density of the 
heaviest L-dyons. Their results are plotted in 
Figj8] In a more recent paper by Bornyakov 
and co-workers [30] an improved identification 
of the dyon type was preformed by adopting 
the fermion method and also by computing the 
value of the Polyakov loop at the dyon center: 
it should be 1 for the M dyons, but —1 for "ad- 
ditionally twisted" L dyons. 



19 



0.6 -i 




FIG. 8: (color online) The density of the topological clusters (in fm~ 4 ) versus the temperature T/T c 
of the SU(2) pure gauge theory. The open blue circles show static dyons, identified as M-type, while the 
open red diamonds are for calorons or L-type dyons [29]. The closed diamond at T/T c — 1.5 is from [30] . 
in which L-type dyons were identified directly by fermionic zero-modes and the value of the Polyakov loop 
at its center. The dashed and solid lines correspond to semiclassical expectations for M, L dyon density, 
with parameters defined in the text. 



Their result is shown in the same figure by 
the filled diamond, and its suggests the M 
density about twice larger than in the other 
method. Unfortunately, the latest study was 
done only at one temperature T = 1.5T C . The 
factor two discrepancy between the results in 
Ref. [29 and in Ref. [30 is due to different ef- 
ficiency of their topology-searching algorithms. 
The values of the the paper Ref. [29] mean to 
be the lower bound on the densities, as the au- 
thors also report also rather large number of 
"unidentified topological clusters" in the same 
table. 

It is instructive to compare the dyon densi- 
ties observed on the lattice with the analytic 
semi classical expressions. Those include the 
holonomy v(T) = v(T)/(2ttT) and the effective 
action So = 87r 2 /g 2 in a combination 

n M = CV 8l//3 exp(-S ^), (D5) 

where C is the normalization constant and 
the product of the pre-exponent powers of the 
holonomy is taken from the DGPS weight for 
the LM pair [24 . The density of L + L-dyons 
is given by obvious change v — » v = 1 — v. Since 
the holonomy v(T) goes to small values at high 



T more rapidly than the charge is running, one 
expects a significant excess of the M-dyons over 
the L-type ones at high T. The difference be- 
tween the two densities should however disap- 
pear at T = T c , with a characteristic singularity 
related to critical behavior of the SU(2) theory. 

As it is clear from FigjH], the observed depen- 
dence shown by curves can be well reproduced 
by the semiclassical expression, with two fitted 
parameters: 

C = 1.6/m" 4 ,5 =3 (D6) 

We also note that the point from [30] at T = 
3/2Tc w 450 "MeV" corresponds to rather 
small dimensionlless density 

n L+ i/T A « 0.007. (D7) 

or volumes/dyon above even the highest values 
used in our simulations. 

We recall that the lattice data apply to 
a fully interacting ensemble, e.g. including 
the screening part of the action and that the 
data/expressions discussed so far apply only to 
a pure gauge theory, i.e. to Nf = 0, which 
is characterized by the weakest coupling along 
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the critical line. Hence, by increasing num- 
ber of quark flavors we expect to assess a sig- 
nificantly stronger coupling domain, hence a 
smaller So ~ 1. If so, eventually there would 
be no numerical suppression and the diluteness 
should become n LJr i/T A ~ 0(1). Such a dense 
liquid of (non-semiclassical) dyons is obviously 
not expected to be described by the simple 
formulas used above. This regime should be 
the subject of future theoretical and numerical 
studies. 

While grouping of the dyons into instantons 
suggests about equal number of L, M dyons, 
this equality is only reached at and below T c , 
while at T > T c the lighter M dyons are more 
numerous than heavier L ones. The relevent 
lattice data are shown in Figj5J As one can see, 
the density of M-dyons grows at high temper- 
ature, as their action decreases. The density of 
L-type, associated with calorons, decreases, at 
large T. We emphasize again that the identified 
topological objects in the lattice simulations are 
just a sub-set of the of the observed topological 
clusters. 

Only the solid diamond corresponding to the 
analysis of Bornyakov et al [30] should be con- 



sidered as a quantitative determination of the 
L-type dyon density. Indeed, this was done by 
combining the fermionic quasizero-modes filter- 
ing procedure with the action smearing to lo- 
cate topology, and by computing the sign of 
the Polyakov line at the center to identify the 
L and M— type dyons. Such a work focused 
on T = 1.5 T c SU(2) gauge configurations in 
a 20 3 x 4 lattice. These authors confirmed 
that several main predictions of the semiclas- 
sical theory remain valid in their lattice config- 
urations. In particular, all fermionic zero and 
near-zero-modes are locally chiral ( a general 
property of their topological nature), the peri- 
odic fermions interact with M and antiperiodic 
fermions with L dyons, as expected from semi- 
classical solutions. The L dyons are dilute and 
paired with L ones into clusters" discussed in 
[22] . The chiral symmetry of the antiperiodic 
fermions is unbroken. The M-type dyon en- 
semble is dense and thus chiral symmetry for 
periodic fermions is broken. They also found 
that while the L-type dyons occupy only about 
3% of the lattice volume, the M dyons occupy 
a significant fraction of it. 
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